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ABSTRACT 


Theoretical implications of the design of a two-level maser of the 
molecular-beam type are presented. The design and operation of a 
microwave maser that employs a beam of ammonia molecules are 
described. 


I. INTRODUCTION 


A NOTABLE advance in the technique of molecular-beam spectroscopy in 
recent years has been the development of a device called the ‘‘maser” 
(a device involving a beam of molecules which gives ‘“‘Microwave Amplifica- 
tion by Stimulated Emission of Radiation’’) by Gordon, Zeiger and Townes,! 
who employed a molecular beam of ammonia in their first experiment. 
The gas maser is the precursor of the solid-state masers?»? in which para- 
magnetic crystals are used. These masers are already being put to use in 
radio astronomy. Because of a fairly good collimation of the beam, the 
Doppler broadening is small; the emission line is approximately 7 kc. wide. 
The mean frequency deviation of an oscillating maser is quoted as a few 
parts in 10%. The gas maser provides a frequency standard of high preci- 
sion and thus it can be used for operation as an “atomic clock”. Aside 
from its use as a spectrometer, the gas maser has recently been employed 
to secure a direct experimental verificationt of one postulate of the Special 
Theory of Relativity, the constancy of the velocity of light. 


* This work, which is based on an Sc.D. Thesis submitted by H. G. Venkates to the Depart- 
ment of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts, U.S.A., 
August 25, 1958, was supported in part by the Ministry of Education, Government of India, and 
the Indian Institute of Technology, Kharagpur, India. 


¢ Present Address: Department of Physics, Indian Institute of Technology, Kharagpur, 
India. 
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The principles of operation of this device can be extended by using 
other molecules, for example, alkali halides. There must be some modi- 
fication in the construction of the beam sources, but the range of frequency 
standards® will thereby be extended. It is also conceivable that the device 
may provide a means of verifying a prediction of the General Theory of 
Relativity concerning the influence of the gravitational field on time. Two 
spectral lines generated by two identical atoms (or molecules) in the same 
quantum trasition, one atom being situated in very high gravitational 
potential (for example, close to the sun), and the other in lower gravitational 
potential (for example, on earth), differ in frequency. The radiation reach- 
ing the earth from the distant atom has a longer wavelength than the radia- 
tion generated on earth. , 


The operation of the maser relies on the possibility of obtaining a 
beam of molecules, virtually all of which are in a quantum state of energy 
higher than other available states and from which they may be induced to 
make a transition to a state of lower energy. The balance of the energy 
appears as radiation. The “‘state selection” of the molecules has been 
accomplished in the maser by means of an inhomogeneous electric field that 
sorts the molecules in the beam, according to their quantum states, as they 
traverse it. The state-selected beam then enters a tuned microwave cavity 
in which, under the influence of a radiation field in the cavity, the beam 
emits power. If a state-selected beam of sufficient strength is sent into 
the cavity, the radiation emitted by the beam will be of sufficient intensity 
to compensate for the loss in the cavity wall, and a stationary wave pattern 
can be set up in the cavity. The cavity is then said to be functioning as 
an oscillator. 


The state selection by the static fields depends on the differential Stark 
effects on the molecules in the several quantum states. The r-component 
of the force exerted upon a molecule is given by 


where W is the potential energy of the molecule in the electric field E. The 
force F, on the molecule depends, therefore, not only on the field gradient 
but also on the quantum state of the molecule (which determines W for a 
given electric field). Hence, the separation of molecules is brought about 
in the inhomogeneous field. A molecular beam of ammonia having iso- 
lated pairs of rotation-inversion states (each pair of levels is described by 
the rotational quantum numbers J and K, with the inversion index sup- 
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pressed) is ideal in this respect, since the forces associated with the two 
states are opposite in sign for limited values of the field. Thus 


p?M?K2 

J?(J + 1)? dE (2) 

(yt ivseee| 
2 Y y+ IP 

where j is the electric dipole moment of the molecule; vp, the inversion 
frequency; J, the quantum number for the total angular momentum; and 
K, M, the quantum numbers obtained from the projection of J along the 
molecular axis, and along a space-fixed axis. If the field structure has axial 
symmetry, and the field increases with lateral displacement from the axis 
up to a value E,,,,, then the molecules in the inversion state of higher energy 
emerging from a source that has small angular divergence about the axis 
of the field are subject to a lateral force inwards and can be trapped into 
a microwave cavity (Fig. 1), provided that their transverse velocity does 
not exceed a critical value v, given by 





F, = F 


We = 4m’, 


where W, is the Stark energy at the surface of the focuser. The emitted 
power and the (angular) frequency of oscillation w, of the stationary waves 
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Fic. 1. Diagram showing the essential parts of the maser oscillator. 


in the cavity, slightly modified by the lack of perfect tuning of the cavity 
to the resonant (angular) frequency w, of a particular inversion transition, 
are given.” by 


2 (N’ L\* l 2 sin? 
Pp = 4G? (N’tw) (7) (1 + loz a7) (®) Ge" (3) 
@— W __ We — %,2Q 1 
Wo = We Wo 7,” (4) 


where 


1 1\3 
6= FATAC + log sr) [(w; — wo)? + G?}} ; (5) 
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7, is very nearly equal to R/ve (see Appendix); N’ represents the flux of 
molecules in the cavity; G, the electric field amplitude in the cavity; J, L, 
the lengths of the focuser and cavity, respectively; R, the cross-section 
radius of the focuser and cavity; we, the resonant frequency of the cavity; 
and Q, the quality factor of the cavity. The oscillation condition is given 
by 


1 4y/m Np? f< 

Of & a ve@L (1 + log sr), (6) 
where N represents the excess of molecules in the inversion state of higher 
energy in the incident beam; and a, the most probable velocity of the mole- 
cules in the beam source. 


II. GENERAL EXPERIMENTAL DESIGN 
To find the operating characteristics of the maser explicitly, we write 
the critical velocity occurring in the oscillation condition as follows: 
‘ M?K2 4 
We = 4 mv,.?2 = [ (ave)? + 42%... »* POL 7] — thy. (7) 
The oscillation condition can then be written 
E> Gt DG 1, tm, 
_ p2mt (QA)? N# 2hvp?/m QA N 
(8) 
where 


32m)! l 
A= Span L(1 + loe sr). ' 


The terms in the braces in (8) decrease rather rapidly for N > 10". 


Rather than the maximum field of the focuser, the potential given to 
the rods of the focuser is actually measured. The relation between the 
focusing voltage V and the field E is determined by the geometry of the 
focuser. In the actual experimental arrangement a 2p-pole focuser® was 
used. This consists of 2p-parallel rods of metal arranged symmetrically 
on a cylindrical surface. The field arising from the geometry of the focuser 
has a cross-section component given by 


CpV (5) 


[1 +(5)" —2(¢) cosas] 


where r is the radial displacement, and ¢ is the azimuthal angle. 


E(r, ¢) = 





(10) 
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Furthermore, 
a oo Po np 
2-8) 
Cc Sp — Qn ~ 1) ’ 


. a ss * 


b= + [0 (& — po], 


where p, is the cross-section radius of the rods, & that of the focuser and 
oo <&. The relation between the maximum field and the voltage is given 
approximately by 


oe 
| iia 2&, V, (11) 


with the condition py < & that was assumed in the derivation of (10). Com- 
bining (11) with (8), we can write the following relation for the minimal 
oscillation condition, which is valid within a limited range of N: 


V, 2£5 JG + 1) (hyo) aS Tt 
o™ p * MK * gmt (QA)? N? 
1 1 7 
x (1+ gaye QA: Wt)» (12) 





where V, is the starting voltage for producing self-sustained oscillations in 
the cavity for a given density N (in the beam source) of the molecules in 
the inversion state of higher energy. The oscillation condition depends, 
therefore, on the cross-section radius of the focuser, the lengths of the cavity 
and focuser, the properties of the molecule (m, ji, v9), the quantum state 
of the molecule (J, K, M) and the thermal conditions of the source (N, T). 
The estimated values of V, for two values of the density are: 


‘N =10"; 104 
V,=9 kv.; 2-15 kv. 


We shall now describe the cavity. There are two properties of the 
cavity to be considered: the quality factor Q, and the mode of the electro- 
magnetic oscillation in the cavity. If cavity tuning is to be accomplished 
by varying its diameter alone, then the TMo;>) mode is the appropriate one. 
This mode is also associated with a high figure of merit. The frequency of 
the mode is determined as v = 2-405 c/7D, where D is the diameter of the 
cavity, and cis the velocity of light. The transverse and longitudinal com- 
ponents of the field in a circular wave-guide are indicated in Fig. 2 (a). 
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Figure 2 (5) illustrates the coupling of the cavity to an external wave-guide 
system through an iris located halfway between the ends of the cavity. The ex. 
ternal Q,, which measures the loss of power through the iris, and the internal 
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(b) 
Fic. 2. (a) Field configuration in a circular wave-guide. (6) Coupling of wave-guide section 
to cavity. 
Q,, which measures the dissipative loss in the wall of the cavity, are both 
experimentally determined (with low precision) as 


: a ae Ss ea l 
S“rrioie* Re ee | Q-Q* 0,’ 


where v is the centre frequency; 4v, the width of the cavity mode reso- 
nance; and r, the reflection coefficient. It is desirable to have the external 
Q, large compared with Qo, in order to produce oscillations in the cavity 


most easily. 





Ill. DETECTION SYSTEM 


The method of detection depends on the amount of output power. 
This can be estimated,® under saturation conditions, from the formula 
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where W is the energy stored in the cavity. With the typical values A = 
1om.2, /= 20cm., L = 10cm., ve = 10° cm./sec., Qe = 100,000, we obtain 


Pin. ¥ 3 X 107° watt. 


The noise power that is in competition with the signal power from the cavity 
is KTb, or 10-1 watt, with b equal to a few megacycles. The observation 
of the signal is feasible if we employ a superheterodyne method of detection. 
A block diagram of the microwave circuit used in the actual experiment 
is shown in Fig. 3. The power from the signal generator, tuned to the 3-3 
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Fic. 3. Block diagram of the microwave circuit. 


inversion frequency, is transmitted to the absorption cell containing ammonia 
gas at very low pressure (50 microns) and to the discriminator, by using 
a magic Tee. A directional coupler is used to transmit part of the power 
to the maser cavity through the attenuators. The power reflected from 
the cavity is transmitted to a balanced mixer. The output from the local 
oscillator (a second klystron tuned 30 mc. away from the frequency of the 
signal generator) is transmitted to the balanced mixer, and also to the dis- 
criminator. These two signals are mixed in the rectifying crystals of the 
balanced mixer and a signal with the 30 mc. frequency difference is amplified 
by the intermediate-frequency amplifier (IF strip). The output of the IF 
strip is exhibited on the oscilloscope. The signal generator is modulated 
and the local oscillator stabilised. The discriminator in the circuit is de- 
signed!® for the purpose of stabilising the local oscillator or the signal gene- 
rator. 
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IV. DESIGN OF AN AMMONIA MASER 


Figure 4 illustrates an ammonia manifold that is the source of the mole- 
cular beam. Figure 5 shows the parts of the maser assembly in the vacuum 
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Fic. 4. Gas manifold. 


chamber. The manifold consists essentially of an ammonia tank, a needle 
valve, and a bubble chamber. The system is evacuated to a pressure of 
10-3 mm. by using two forepumps. The stop-cocks E, G are provided to 
isolate the manifold from the pumps in operating condition. The primary 
tank contains ammonia at very high pressures (several atmospheres). This 
pressure is scaled down to one atmosphere in the secondary tank. The 
ammonia gas bubbles through the Octoil in the chamber when the stop-cock 
D is closed and the needle valve C is slightly opened. Octoil ‘‘S” is used 
in the bubble chamber because it has a low vapour pressure of 10-* mm. 
at room temperature. 


The main vacuum chamber, 30 in. long, and Sin. in diameter, is made 
of stainless steel. It is closed on top with a lid that supports the main parts 
of the maser assembly, which are held in position by steel mounts with 
Teflon spacers. To prevent the vapours of the oil in the diffusion pump 
at the lower end of the chamber from condensing on the parts, a baffle is 
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provided. The whole chamber is provided with a liquid-air jacket that is 
welded to it. The jacket is an evacuated double wall; it is used to cool 
the internal parts of the chamber, especially the fin-like condensers in thermal 
contact with the inside wall. 


The gas from the manifold enters the beam source in the chamber. 
The beam source is a cylinder of stainless steel, 3 in. long and 0-5Sin. in 
diameter. It has a nozzle that is constructed as follows: a long ribbon of 
nickel, 0-015 in. wide, is fed between interlocking geared wheels so that 
the ribbon is transversely corrugated. The corrugations are parallel grooves 
of semicircular cross-section shape, with a diameter of approximately 5 mils. 
When alternately corrugated and uncorrugated strips are stacked, narrow 
channels are formed between the layers in a rectangular array that has a 
honeycomb structure to serve as a nozzle. Since the beam source is likely 
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to cool to a very low temperature when the vacuum chamber is cooled to 
liquid-nitrogen temperature, provision must be made for its heating in order 
to prevent condensation of ammonia. The beam source is an independent 
unit that can easily be taken out of the chamber and replaced by another 
of suitable design for work with other molecules. 


The focuser consists of 8 rods of stainless steel, approximately 0-062 in. 
in diameter. The rods are spaced equally on a cylindrical surface and 
supported by Teflon spacers. The cross-section radius of the multipole 
focuser is 0-2 in. and its length is 10in. The rods are sufficiently far apart 
to support a potential difference of 30kv., or so, without sparking. The 
two sets of alternate rods are connected to thick copper wires that are taken 
out of the vacuum chamber through glass insulation inside the supporting 
steel tubes on which the maser parts are mounted. 


The cavity is designed for the TMoi) mode. It is a cylindrical cavity, 
open at both ends, mide of coin silver. The cross-section diameter is very 
close to 0-378 in. at room temperature, and it is 4-74in. long. The dia- 
meters of the cavity and iris are brought to the desired dimensions by 
observing resonant frequency and external Q, and correcting them by lapping 
the cavity. The wave-guide connecting the cavity to the external system 
through the iris has its shorter side parallel to the axis of the cavity. It is 
closed by a mica witidow so as to be vacuum-tight. Exact tuning of the 
cavity can be accomplished under actual operating conditions by electrical 
heating. For this purpose, the cavity is enclosed firmly in a lavite block, 
on which a bifilar heater of nichrome wire is wound. Bifilar winding is 
essential to exclude from the cavity the magnetic field that will affect the 
resonant frequency of the molecule. The focuser is surrounded by the 
structure shown in Fig. 6. This consists of a set of copper discs, each having 
a circular opening, approximately 2in. in diameter, through which the 
focuser passes. The purpose of introducing the baffles is to provide a large 
cold surface on which extraneous molecules, especially the inactive ones 
scattered by the focuser, can be trapped. 


The alignment of the beam source, the focuser, and the cavity is optically 
tested to ensure that the molecules emerging from the beam source actually 
pass through the focuser and the cavity. 


VY. MASER OPERATION 


During experimentation we should aim at achieving a pressure of 
10-* mm. in the vacuum chamber when the molecular beam is on. The gas 
manifold is also evacuated to a pressure of 10-°*mm. The stop-cock F is 
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open to the main chamber, and G, E closed. Then the control valve is 
closed, and the gas from the primary tank admitted into the secondary tank 
at a pressure of approximately one atmosphere, as indicated by the mano- 
meter. The ammonia gas can be made free of any possible contamination 
as follows. The stop-cock B is closed, and the secondary tank cooled to 
liquid-nitrogen temperature. The gas crystallises in the tank. Now, C 
may be opened. The presence of residual gases in the tank is indicated by 
a slight rise in the pressure, as indicated by the ion-gauge meter. These 
gases are pumped out, and thereby the ammonia gas is purified. The valve 
C is closed again, and the crystallised ammonia exposed to room temperature. 
When the stop-cock D is closed and the control valve slightly opened, the 
gas bubbles through the Octoil at a uniform rate. The flow rate determines 
the number of molecules per second that effuses into the vacuum chamber. 
With tubing of 0-25in. diameter in the bubble chamber, one bubble per 
second is the optimum rate; this rate can be varied conveniently during 
operation. Before the gas is allowed to enter the vacuum chamber, the 
chamber is cooled to liquid-nitrogen temperature, and then the cavity is 
heated so that it resonates with the 3-3 inversion line of ammonia. 


Cavity tuning is accomplished as follows. The power from the signal 
generator, transmitted through the absorption cell, is monitored. The 
wavemeter is set to the frequency of the 3-3 line (23,870 mc.), which is 
correct within a few megacycles. (This approximate setting is made to 
ensure that we are in the proper region.) When the klystron is tuned to 
this frequency, the klystron mode shows the pip of the wavemeter as it is 
observed on the oscilloscope. When ammonia is introduced into the cell 
at a pressure of approximately 50 microns a sharp absorption line appears. 
The wavemeter is set to the absorption line, as well as it can be, visually. 
At this stage, the y-input to the oscilloscope is taken from one of the crystals 
of the balanced mixer so that the signal is the klystron mode (reflected from 
the maser cavity) with the wavemeter pip. For this observation it is neces- 
sary to use a 20-db. transistor amplifier. The resonant frequency of the 
cavity is made to coincide with the wavemeter setting by carefully adjusting 
the heating current; thereafter the wavemeter is detuned. It will be neces- 
sary to tune the cavity slightly around this setting when an actual search 
for the emission line is being carried out. 


Now the IF strip is brought into the circuit. The output of the local 
oscillator at a frequency of 30 mc. away from inversion frequency is trans- 
mitted to the balanced mixer. The output of the IF amplifier now becomes 
the y-input of the oscilloscope. The local oscillator is stabilised at a fre- 
quency v)»+ 30mc., as follows. The output of the discriminator crystals 
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is monitored. The discriminator cavity is tuned until the cavity-mode signal 
appears on the oscilloscope. The symmetrical discriminator pattern is 
obtained by balancing the outputs of the two discriminator crystals. The 
sweep is then removed and the local oscillator stabilised. 


The cavity is now set in the narrow region about the inversion frequency 
vo, the local oscillator is stabilised, and the signal generator is swept about 
the inversion frequency. The output of the IF amplifier is seen on the 
oscilloscope as a spike. By slightly tuning the discriminator cavity in the 
local-oscillator circuit, the spike is observed to pass through a minimum 
that indicates the setting of the maser cavity. The power transmitted from 
the signal generator to the maser cavity is attenuated approximately 80 db. 
The control valve is then opened slightly and the gas allowed to bubble 
through the Octoil at a rate of one bubble per second. Voltages of from 
+ 5kv. to + 15kv. are applied to the focuser rods. The sweep on the 
signal generator is gradually cut down while the signal is kept in view. 
Finally, the signal generator is stabilised, and the local oscillator swept to 
exhibit the maser oscillation signal on the oscilloscope. 
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thank the Ministry of Education, Government of India, Dr. S. R. Sen 


Gupta, Director of the Indian Institute of Technology, Kharagpur and the 
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Institute of Technology. 


REFERENCES 


1. Gordon, J. P., Zeiger, H. J. Phys. Rev., 1955, 99, 1264. 
and Townes, C. H. 


Scovil, H. E. D., Feher,G.  JIbid., 1957, 105, 762. 
and Seidel, H. 


Strandberg, M. W. P. Ibid., 1958, 109, 1988. 
Davis, C. F. Jr., Faughnan, 
B. W., Kyhl, R. L. 
and Wolga, G. J. 


Cedarholm, J. P., Bland, Ibid., Letters, 1958, 1, 342. 
G. F., Havens, B. L. 
and Townes, C. H. 


Peter, M. and Strandberg, “Theoretical and experimental study of molecular-beam 
M. W. P. microwave spectroscopy,” Technical Report 336, Research 
Laboratory of Electronics, M.I.T., August 26, 1957, 





H. G. Venkates and Proc. Ind. Acad. Sci., A, Vol. LI, Pl. I 
M. W. P. Strandberg 


ee gee ee 


Fic. 6. Photograph of maser assembly. 























Design and Operation of a Molecular Oscillator 135 


Venkates, H. G. and 
Strandberg, M. W. P. 


Shimoda, K. 


Townes, C. H. and 
Schawlow, A. L. 


Vonbun, F. O. 


Wolga, G. J. 


“Operating characteristics of a molecular-beam maser” 
(to be published in J. Appl. Phys., U.S.A.). 


J. Phys. Soc. Japan, 1957, 12, 1006. 


Microwave Spectroscopy, McGraw-Hill Book Company, Inc., 
New York, 1955. 


‘*A multipole separator for masers,’’ Engineering Report, 
NRE-1203, U.S. Army Signal Engineering Laboratories, 
Fort Monmouth, N.J., April 1957. 


““Microwave resonance studies of absorption and hyperfine 
interactions of F centers in alkali halide crystals,” Ph.D. 
Thesis, Department of Physics, M.LT., August 1957. 








136 H. G. VENKATES AND M. W. P. STRANDBERG 





APPENDIX 
The time 7,, introduced in Eq. (4), for the frequency pulling is given® 
by 
f f yer dx dy +f f sae dx dy 
™1= ; . y=0 sny vy=0 s=0 - 
f f Je dxdy +f f pom dx dy 
y=0 s=y v=0 s=0 
where 
Jit, »« £ 
o-E5* °“gro: 


The integral sum in the numerator is evaluated as 


o(-)+ 3 alata ove — Lam 
~0-8)-Gie-he-49) 


+ &(5 4% — ge—5at)+... 


Since, in practice, a« 1 and b« 1, we may neglect the third and higher powers 
and write for the integral sum in the numerator 


501-5) an e('- 


The integral susa in the denominator has already been evaluated.® It 
is equal to 


2R(x) (1 + ler). 
The time 7, is then given by the expression 
,;—2L 
oe 
C1+ log + 


In. the experimental design / = 2L, and hence we may write 
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INTRODUCTION 


1. Notations—The purpose of the present paper is to show that some 
results in the theory of geometrical difference equations can be proved by 
the use of the matrix calculus. In particular, the so-called regular case may 
be treated completely whereas the irregular case brings some difficulties 
which cannot be overcome. It may be mentioned that the analytical be- 
haviour of the solutions of a geometric difference equation has been treated 
completely by Trjitzinsky in a very extensive paper.! 


Let g be a real parameter, 0 << q< 1. The operator of the geometrical 
difference is defined by 


f(x) LL); or f = 0 (6 f)(r = 1, 2, ...). (1.1) 


The inverse operator may be expressed by means of an infinite series. If 
OF (x) = f(x), we get? 


F (x) — F) =(1—9)x E gi fxd). (1.2) 
The right-hand side is sometimes written as 
S f(x) d@, »). (1.3) 


The general linear geometrical difference equation may be given in the form 
by = A(x) y+ z, (1.4) 
y (x) being a column vector, z being a given vector and A (x) a given matrix. 





1 Cf. (1). * Cf. (2). 
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The equation 


Ag (x) O f + ay (x) POF + «6. + Ana (x) OF + an (XS = 8 (x) 
(1.5) 


may be regarded as a special case of (1.4), if we introduce the vector y with 
the components 
W=f, Yo=A,.--5 n= Of 
the vector z with the components 
Ze... = Zy, = 0, =F 


43> 


and the matrix 


0 l 0 0 
[° 0 l 0 
A= . ee ; 
0 0 Pee er ae oe l 
4m 41 Gra | 

Ay Ay ay 


Because of (1.1), the equation (1.2) may be brought into the form 


y (qx) = B(x) y(x) +2 (1.6) 


where 


B=I—x(l —q)A,z = — x(1 — q)z (1 the unit matrix). 
Let the function eg (¢) be defined® by 


Re aed, (1.7) 
IT (1 — q't) 


r=0 





This meromorphic function may be represented by the power series 


ee ee aoe See Te 
‘I—q' @—-q:' °° "d—an 





ok. (1.8) 


where 


(1 — qn = (1 — g) (1 — g’)... A — q”). (1.9) 
The series (1.8) is convergent for |¢| <1. 


8 Cf. (2). 
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Let P(t) be a matrix, bounded for |¢|<%. Then the matrix 

F [P (¢), t] = (I—¢P (0)? (I—@tP (qn [F — @2tP (g’)}" ... (1.10) 
has a meaning for sufficiently small |¢|. This may be proved as follows: 
Surely, a certain | t, |<) will exist so that for | ¢|< t, each factor 


I — g™tP (q"t) (7 =0, 1, ...) 
is not singular and that its inverse may be written as 
I+ g™tQ (g""), 


Q(t) being bounded too. Let m be a bound of the elements of Q and let M 
be the matrix with all its elements equal to m. Then each element of the 
matrix (1.10) will be majorized by the elements of the matrix 

(I + *M) (I + gtM) (I + gM)... (1.11) 


Now, this matrix may be developed into the everywhere convergent power 
series* 





tM 1M reraaya_ 2M" 
itpagttiagt tea gt OD 


and, therefore, each element of (1.10) has a convergent majorant. That 
means the existence of (1.10) has been proved. If P is a constant matrix 
all the factors of F (P, ) commute inter each other. Then it is possible to 
develop F (P, ¢) into the series 
¢*ps 

:*c a 13) 
which is convergent for sufficiently small | ¢| and with regard to (1.8) may 
be denoted symbolically by eg (tP). 





14,2 


2. The Regular Case-—We regard the homogeneous equation (1.4), 
1.€., 
dy = A(x) y (2.1) 
and assume that y (0) 0 exists. Then, A(x) must be bounded for x = 0, 
otherwise the equation [cf. (1.6)] 


y (gx) =(I— x — gg) AG@)] y@) (2.2) 
had no meaning for x = 0. Now, one can see immediately that the vector 
y = F[A(), x(1 — g)] yO) (2.3) 


is a solution of (2.2) and also of (2.1). Consequently, the general solu- 
tion of (2.1) is® 





‘© Cf. (). * Cf. (I). 
Ag 
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F[A (x), x(l—g] R(x) y), (2.4) 
where R(x) is q-periodical, 
R (qx) = R (x). 


The solution (2.3) is regular in a certain neighbourhood of x = 0. 


There exists one more way to find the solution of (2.1). We consider 
the matrix 


S=1+(1—g)x EF iA (qix) + (1 — 9)*x* F QA (xq/) 
j=uo $20 
x S gkA (xqit®) + (1 — g)*x3 5 g*A (xqi) 
k=0 {30 
x F gta (xqit®) 3 gta (xgit*t +... (2.5) 
k=0 tso 


If A is bounded, the convergence of this series can be proved in the 
same way as we proved the convergence of (1.10). Replacing A by a 
matrix M with all its elements equal to m in (2.5), we get 


[+(1—q)xMZYqi + (1 —q)? x? M?Z q4 F gk 
+ (l=) x? M9 Eq) Tq Faq +... 
= €q [x (1 — q) MJ 
which is therefore a majorant for (2.5). 


Now, we can see immediately that 


6S = A(x)S (2:6) 
and consequently 
¥ (x) = Sy 0) (2.7) 


must be a solution of (2.1). Because $(0)=I, the solutions (2.7) 
and (2.3) must be the same, and we see in the same time that 


This formula may be interpreted as an expansion of the infinite matrix 
product (1.10) which is analogous to the well-known 2-matrix arising from 
A(x) by repeated integration. 


4 


3. The Inhomogeneous Equation —To find the general solution of the 
equation 


Gy=Al(x)yr2z (3.4) 








——ae 
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we need a particular solution y+ and then we get the general solution add- 
ing this y+ to the general solution (2.4) of (2.1). To get the particular 
solution, we apply the method of the variation of parameters. We set 


y* = F[A(), x(1 — g)] w@) (3.2) 
where w(x) is another unknown vector. Then we get 


dy* = A(x) F[A (x), x (1 — g)] w(x) + F[A x), x9 (1 — g)] Ow. 
(3.3) 


Substituting (3.2) and (3.3) into (3.1) we get 
z= F[A (qx), gx (q — 1)] Ow (3.4) 


and we can find w by means of the formula (1.2) or (1.3) as 
w= SF[A (qx), gx (1 —g)}*z(x) dq, x) (3.5) 


Therefore, the general solution of (3.1) has the form 
y=F[A(a), x(1 — 9)] 


x[ R (x) vy) + S F {A (gx), gx (1 — g)}*z(x)d (@x)| 
(3.6) 


4. Constant Matrices.—If the matrix A in (2.1) is a constant one we 
can represent the solution by means of the function eg (#) in the same way 
as the solution of a differential equation may be composed out of exponen- 
tial functions. 


Let P be a non-singular constant matrix which transforms the matrix 
A into its canonical form J. 


PAP = J. (4.1) 
Let 

poe Oe (4.2) 
then we have 

du = Ju. (4.3) 
Applying (2.3) and considering (1.13) we find } 

u(x) = eqg[x(1 — g) JJu@) i (4.4) 
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and therefore 


v(x) = Peg [x (1 — g) J] Py (0). (4.5) 


Now, J being a canonical matrix, eg (tJ) can be calculated very easily. Let 
T be one of the “boxes” of J of the order k, belonging to the eigenvalue 
a of A, 


aio0 . 60 
Set... o¢ 
000 a 1 
0 0 0a 
Then we have 
P me in a 
a= 28 qr 4.7) 


Calculating the powers of T we find for r = 0, 


ties 


im 3 (er Ger Gem 
0 a’ C es Gi: oe DT ee (4.8) 


0 0 0 ee a’ 


Substituting (4.8) into (4.7) we see that all clements in the main diagonal 
of (4.7) are 


LJ 


>, t* Tr 


r=0 


The elements in the first oblique row parallel to the main diagonal are 


“——  -., 
ar a dines 
Therefore 
a (at) teq'(at) ... (; ul 1) thle gtk-w (at) 


eq (tT) ‘\ 0 eg(at) ... (, 2) tk-2¢,(k-2 (gt) 


: 0 0 di &q (at), 
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Now, @q (tT) is a “box” of eg (tJ) in the same way as T is a “box” of J. 
Consequently, from (4.5) we learn that the general solution will be a linear 
combination of 

eq (x (1 — g) ai], xeq’[x(1 — Qa], ..., xP tegt*e [x (1 — 9) ai] 
where the ay are the eigenvalues of A belonging to different elementary 
divisors and the k; the corresponding orders. We see that the solution has 
an analogous shape as the solution of a linear differential equation, the 
exponential functions having been replaced by eg(t) and its derivatives. 
But because eg (¢) is a meromorphic function the solution of (4.4) will be 
a meromorphic function also whereas the solution of a linear differential 
equation with constant coefficients is an entire function. 


5. The Singular Case.—Now let us weaken the assumption that y (0) 
will be a constant finite vector, different from zero. We are restricting 
ourselves to the homogeneous case. We consider the equation in the more 
general form [cf. (1.6)] 


y (qx) = B(x) y (x) (5.1) 
and assume that B(x) may be bounded and that 
B (0) = C (det C¥ 0). (5.2) 


Without loss of generality we may assume that C is given in its canonical 
form otherwise we may apply a linear transformation like (4.2). Now we 
write 


y (qx) =[C + xD (x)] y () | (5.3) 
and introduce 7+ by 
a _ logx 
x=q* or ”=To8q" (5.4) 
Let the matrix T given by (4.6) be one of the “boxes” of C and let 
a= q* 
The matrix 
x* x#q-* () xg~2a (3) = xagi-ke ( bed 1): 
Na(x) = 10 x xq*(7) ibe xegitke ( bd 2) 
0 0 ie a 


(5.4) 
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satisfies the geometrical difference equation 
Na (gx) = TNa (x) (5.5) 


which may be verified easily. Composing the different “boxes” N, to a 
matrix N of the order n we find that N will satisfy the equation 


N (qx) = CN (x). (5.6) 
Now we set 

y (x) = N(x) u(x). (5.7) 
Then u(x) will satisfy the equation 

u (qx) = (I + xN-!C-?DN) u (x) (5.8) 


which is of the special type considered in 2 if the matrix G = xN-1C""DN 
is bounded. But this condition is not fulfilled at any rate because the 
elements of G contain powers of x multiplied into powers of logx. A 
sufficient condition for boundness will be that all the occurring powers of 
x have an exponent larger than zero. If a,,..., ap are the eigenvalues of 
C and if cj = q%, then the exponents which occur in N will be a, ..., ap, 
and the above-mentioned condition may be written as 


1+ a; —o; >0(i, j=1,...,p) 
or 
—t<qgj<0+}. (5.9) 
If it is fulfilled the equation (5.8) can be resolved, u(0) being finite. 


Considering (5.7) we see that the solution y(x) will be weakly singular 
in the neighbourhood of the origin. 


If det C = 0 the equation (5.8) has no meaning, and the discussed 
method cannot be applied. In this case, some of the components of y (x) 
will be strongly singular at x = 0. 


6. Exceptional Points.—In the above reflections, the point x = 0 was 
treated as an exceptional point in which the solutions may get singularities. 
Now, every geometrical difference equation has at least one more excep- 
tional point and that is x = oo. Substituting ¢ = x-! we get an equation 
for y (€) = n (€) which has the origin as its exceptional point, and so we 
have the former case. Sometimes it happens that x = 0 will be a strong 
singularity whereas x = co will be a weak one. 
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In the regular case we can solve the equation (2.1) also for given initial 
values y(a), a0, provided that a is situated within the circle of con- 
vergence of (2.4). In fact, we have 


y (a) =F [AQ@), a(l—g]y@ (6.1) 
and by eliminating y (0) from (2.3) and (6.1) we find 
v(x) = F[A(x), x(1 — q)] [F{A@, a(l — gy @). (6.2) 
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THe C = O stretching frequency in the case of normal ketones, aldehydes 
and esters, occurs at fairly constant position both in Raman effect and 
infra-red.1)2 The carbonyl frequency is determined by the nature of its 
immediate environment. The infra-red studies by Cross and Rolfe* on some 
ketones show that there is an increase in frequency of the C = O line with 
non-polar solvents and a decrease with polar solvents. Their results show 
that the factors which influenee the frequency of the carbonyl absorption 
also influence its characteristic intensity. 

Jesef Goubeau and Ruth Deerdt,* from the study of the influence of 
solvents on the Raman lines of some ketones, have shown that there is a change 
in frequency of all the lines, especially the C = O line. There is an increase 
in frequency with non-polar solvents and polar halides, while the other 
solvents have produced a decrease. 


The present investigation was undertaken to study the influence of various 
solvents on the frequencies and the intensities of the Raman lines in various 
unconjugated and conjugated substances. 


EXPERIMENTAL 


Extra pure samples of the substances distilled at constant boiling points 
have been used in the present investigation. A Hilger E612 two-prism glass 
spectrograph, in conjunction with a Hilger Raman Source Unit was used. 
The spectra of the various substances have been recorded under identical 
conditions. A set of graded intensity marks were also given on the same 
plate and the intensities were computed using a Moll microphotometer. 


RESULTS 


The mixtures used are 5C% by volume. The intensity values given in 
the accompanying Tables I to IV are the ratios of the integrated intensities 
of the various lines to the integrated intensity of 458cm.- line of CCl,. 
Table V contains the values of C = O frequency in various substances in 
both pure and solution states. 
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TABLE I 


C =O Frequency 





x a 
Negus Pure CS, CCl, C,H; CHCl, C,H, Dioxan 
* 
Sub- \ 





stance 

Methylethyl 0-064 - ea 0-054 - 0-035 
ketone 

Ethyl acrylate iy ae a 0-073 es 0-028 


Benzaldehyde 4-540 9-900 0-669 0-344 1-242 0-739 1-820 
Acetophenone 1-153 2+380 0-445 0-175 0°596 0-260 





Benzophenone - 2-390 0-360 - 0-779 0-277 0-786 
Phenyl 0-169 0-208 0-092 0-034 ed a 0-021 
acetate 

TABLE I] 


C =C Frequency 








\ 
oe Pure CS, CCl, CHCl, Dioxan 
Sub- 
stance - \ 
Benzaldehyde 3-472 9-615 0-625 1-477 1-095 
Acetophenone 1-071 3-674 0487 0-679 oe 
Benzophenone -" 4-019 0°474 0-795 1-226 


Phenyl acetate 0°215 0-436 0-088 a oe 
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TABLE III 


C — Ayrom Frequency 

















‘ 

% 

\ Solvent Pure cs, CCl, CHCl, Dioxan 

Sub- \, 

stance \ 
Benzaldehyde 1-520 1-449 0-489 0-601 
Benzophenone 0-973 0-269 0-455 0-530 
Acetophenone 0-407 1-157 0-305 
Phenyl acetate 0-256 1-610 0-311 

TABLE IV 
C — C Frequency 
~~ ? 
\ Solvent Pure CS, CCI, CHCl, Dioxan 
* 
Sub- \ 
stance = \ 
+ 

Benzaldehyde 3-113 3-210 0-364 0-569 0-687 
Acetophenone 0-936 1-633 0-679 0-332 
Benzophenone 1-924 0-215 0-455 1+40 
Pheny! acetate 0-759 1-81] 0-163 0-780 











Ps | 
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TABLE V 


Frequency of C = O Vibration (cm.~) 








x 

*%, Solvent Pure CS, CCl, CHCl, C,H, C,H, Dioxan 
Sub- 

stance \, 

Ethyl acrylate 1720-0 1727-8 1725-0 sia xii 1730-0 


Acetophenone 1680-0 1690-4 1687-0 1680-9 1691-3 1673-3 1681-8 
Benzophenone 1660-0 1665-2 1665-2 1659-7 1664-2 1662-7 1661:9 


Benzaldehyde 1700-0 1704-0 1705-0 1700-8 1708-5 1709-0 1701-7 





DISCUSSION 


The following conclusions can be drawn from the results given in the 
above tables :— 


(1) The frequency of the C = O line shows, in general, an increase in 
the non-polar solvents and no appreciable change in CHC], solution. This 
observation is in agreement with Goubeau and Deerdt so far as the non- 
polar solvents are concerned; but in the case of polar halides they have 
recorded an increase while there is no change in the present investigation. 
The C=C, C —C frequencies do not show appreciable variation. 


(2) The intensity of the C = O line in the unconjugated methyl ethy! 
ketone does not show any appreciable variation in different solvents. 


(3) The intensity of the C = O line in the conjugated aromatic compounds 
is high compared to that in the unconjugated or conjugated aliphatic com- 
pounds even in the solution state. 


(4) The intensity of the C = O line in CS, solution is found to be 
considerably increased whereas in other solvents an appreciable diminution 
is noticed. This diminution is more marked in non-polar solvents. 


In the case of phenyl acetate, where there is conjugation of C = C 
with the benzene ring through the oxygen atom, the intensity of the lines is 
very small compared to the other conjugated aromatic compounds. In 
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this case also the intensity of the C = O line has been considerably increased 
in CS, solution. 


The behaviour of the C= C, C — C and C — H,,,,, vibrations is quite 
similar to the C = O vibration. 


SUMMARY 


The changes in frequency and the intensities of the various unconjugated 
and conjugated substances in different solvents have been investigated. It 
is found that there is an increase in C = O frequency in non-polar solvents, 
The intensity of the C = O line shows in general, consideravle diminution 
n non-polar solvents. The C=C, C —H,,,,, and C—C_ vibrations 
also behave in the same way. 
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A coop deal of work has been done in both Raman effect!~? and infra-red*~ 
on the assignment of characteristic frequencies in carbonyl compounds and 
nitriles. In majority of cases the Raman and infra-red frequencies are very 
nearly the same although occasionally large differences are also observed. 
In the infra-red, some work has been done on the intensity of carbonyl 
absorption*-* of carboxylic acids and also of nitrile absorption®.” but very 
little work has been done in Raman effect. 

In the present investigation, the intensities of the characteristic bonds 
of type C = O, C — H bending, C — Hj;,,, ‘stretching, C —C and C =N 
whose frequencies lie in the region 1700, 1400, 2900, 800 and 2300 cm. 
respectively, have been studied. 


EXPERIMENTAL 


Extra pure samples of the substances distilled at constant boiling-points 
have been used in the present investigation. A Hilger E 612 two-prism glass 
spectrograph in conjunction with a Hilger Raman Source Unit was used. 
The spectra of various substances have been recorded under identical condi- 
tions. A set of graded intensity marks were also given on the same plate 
and the intensities were computed using a Moll Microphotometer. The 
intensity values obtained have been reduced to standard intensity values 
(S) by making use of the expression suggested by Bernstein and Allen.*® 


RESULTS 
TABLE I 
Carboxylic acids 








Nature Standard 
Substance of apes ey " Tse cer. intensity 
vibration _— (S) 
Acetic acid = 1669 0-063 0-15] 
Propionic acid Re C=O 1669 0-073 0-233 
Butyric acid ¥ 1670 0-107 0-434 











TABLE I—(Contd.) 
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Name S‘andard 
Substance _ Of | F —s T/Tese cig inten ity 
vibration , (S) 
Acetic acid 1435 0-124 0-242 
Propionic acid C — H bending 1448 0-207 0-552 
Valeric acid 1430 0-171 0-690 
Acetic acid — 2934 0-211 1-131 
Propionic acid C-H 2956 0-354 2-576 
Butyric acid atigh 2944 0-329 3-009 
Valeric acid 2956 0-484 5-416 
Acetic acid 895 0-325 0-355 
Propionic acid Cc—C 842 0-232 0-316 
Valeric acid 826 0-144 0-295 
TABLE IT 
Carboxylic esters 
Nature Standard 
Substance of Frequency [lass cc, intensity 
vibration — (S) 
Methyl butyrate ‘ 1740 0-104 0-540 
=O 
Ethyl butyrate 1739 0-171 1-025 
Methyl acetate 1455 0-201 0-601 
Amy] acetate - . 1452 0-270 1-471 
Methyl butyrate C—Hbending = 1453 0-230 0-942 
Ethyl butyrate 1455 0-234 1-111 
Methyl acetate 2942 0-163 1-314 
Ethyl acetate 2940 0-183 1-663 
Amy]l acetate C — Ayrios 2940 0-256 3-761 
Methyl butyrate 2940 0-443 4-902 
Ethyl butyrate 2941 0-486 6-242 
Methyl acetate 646 0-161 0-176 
E.hyl acetate c-—C 633 0-125 0-151 
Amy! acetate 620 0-072 0-136 











y 
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TABLE ITI 
Nitriles 
Nature Standard 
Substance of Fr so a T/Tgss cca, intensity 
vibration z (S) 
Aceto nitrile se 2250 0-282 0-852 
Propio nitrile se 2245 0-175 0-750 
CeN 
Butyro nitrile 2244 0-092 0-511 
Benzo nitrile ya 2224 2-119 17-440 
DISCUSSION 


From Tables I and II the following general conclusions regarding the 
characteristic vibrations in both carboxylic acids and esters may be arrived at: 


~ -({1) The intensity of the lines due to C = O vibration is found to increase 
as we go from lower to higher members of the series. 


(2) The intensities of lines in the regions 1400 cm.-! and 2900 cm.-! 
characteristic of C— H bending and C — H,,,, stretching vibrations also 
behave in the same way as in (1). 


(3) The addition of a CH, group in a homologous series has a different 
effect on the low frequency region characteristic of C — C vibration in the 
sense that the intensities show a gradual decrease from lower to higher members 
of the series. 


From Table III it can be seen that (1) the intensity of the line due 
to C= N vibration is found to be decreasing from the lower to the higher 
members of the homologous series. This conclusion is in agreement with 
the observation in the infra-red.’ (2) The large intensity in aromatic 
conjugated substance of benzo nitrile compared to the aliphatic unconju- 
gated nitriles may be explained as due to the combined effect of the general 
increase in intensity to be expected from aliphatic to aromatic and also from 
unconjugated to conjugated substances.*-!° This observation is in conformity 
with that made by Cross and Rolfe® in infra-red absorption. 
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SUMMARY 


The intensities of characteristic vibrations in carboxylic acids, esters 
and nitriles have been studied. In both carboxylic acids and esters the 
intensities of Raman lines due to C=O, C—H _ bending and C—H,,,, 
stretching vibrations are found to increase from lower to higher members 
of the series while there is a decrease in intensity of C — C frequency. The 
intensity of the lines due to C= N vibration is found to be decreasing from 
lower to higher members of the series. The intensity in aromatic conjugated 
nitrile is found to be abnormally great when compared to that in unconju- 
gated aliphatic nitriles. 
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(Communicated by P. L. Bhatnagar, F.N.L, E.A.SC.) 


Let (9;, Pi, 22;) and (gs, p2, 22) be the two flows of an incompressible fluid 
of uniform density p, kinematic viscosity v and kinematic cross-viscosity 
ve, Where g, p, 2 denote the velocity vector, the presssure and the potential 
from which the external forces are assumed to be derived. The two flows 
are said to be superposable, if a pressure (p, + p2 +7) can be found such 
that (q: + G2, Pi + Po +7, 2, + Q,) is also a solution of the Stokes-Navier 
equation with the necessary modification in the initial and boundary con- 
ditions. 


The equations of motion for a non-Newtonian fluid are given by, 


d : ° . 
p(<4 + ui, 545) =ygttt (0.1) 


uy,5 = 0 (0.2) 
1;' being the stress tensor given by 
tt = — pdt + Fydj' + Fadeid;*, (0.3) 
where f* is the body force, 
dj! = uj4 + Wi,5 (0.4) 
is the rate of deformation tensor, F, and F, are functions of some material 


constants and the second and third invariants of the rate of deformation 
tensor d;‘. 


By taking F, = » and F, = pe we recognize F, and F, to be the 
coefficient of viscosity and coefficient of cross-viscosity. 


, 1. For the problem on hand it is convenient to use cylindrical co- 
ordinate system (r, 8, Z) in which z-axis is taken to be the axis of symmetry 
and r is the distance from the axis. The velocity components at any point 
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are given by uy, ug and uz, where the velocity vector q does not necessarily 
lie in the meridian plane, so that 


q = Uriy + Ugly + Uziz, 


where for axial-symmetry the velocity components are functions of r ang 
z only besides the time f. 


(a) Equations of motion —The equations of steady motion in cylindrical 
co-ordinates, with terms of azimuthal variation neglected, are, 


9 
p Uy WU, oa Uu wy at Ug” = dP yr “4. Pre + Pry — Poo (1 1) 
wr '  * (bz r ar dz r , 
_ lu Uy oP oP - 
.: Se : as dz 78 
p (up “6 ~ wo 4. te )= ro 4 4 (1.2) 
p(u yt Mz = ‘> * et yp (1.3) 


along with the equation of continuity 


Uy Ur Wz 

ar ;* =o (1.4) 
where the stress matrix Pj; and the rate of strain matrix e;; in the case of 
non-Newtonian fluids are related as (Braun and Reiner’) 


Prr Pro Prz 1 0 0 Crr &rg &rz 
Por Pog Pez J=— PO 1 OP HHL Cor egg gz 
Pzr Pzg Pzz 00 1 Cor Czy Czz 
Crr rg rz Crr rg rz 
+ He | Cor &¢6 Sez Cor og ez | > (1.5) 
Czr &zg zz Czr Oz Czz 


where ej; = ¢ji 


ou Weg Uy 

wees ~<a 4 
Me dUr Uz 

Cyg = 2 t2= s+ 5 (1.6) 
Uz Ug 

e —_ 2 = e p—— -4 

- dz - 2 


Using the relations (1.4), (1.5) and (1.6), the equations of motion (1.1), 
(1.2) and (1.3) for axially symmetric steady flow become, 
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Expressing the velocity components in terms of the Stokes’ stream 
function %, we have 


1d l dy 
up =, Uy = x (rr, 2), =~; wr’ 
be. 
1) i 
= ; ss ip + Xi, + +(-12 ~ — Plz. (1.7) 


(b) The conditions of integrability.—The conditions of integrability of 
the equations of motion in terms of the Stokes’ stream function % turn out 


to be, 
yp r-*Byb x 
n(t ‘ )—2x 
EN FE 
= vE*Y + 2r¢ 73(% )- 2J 
r 
z 
Lea) _ 1 pay, 1 204 2 (EY 
; EM) r2 dz BY + Fas re 
ay WX X 
rs ae - r) an (= ") jap “ 
Pr OZ 
-3% ent , (1.8) 
and 
l yb rX 
w( 7 
X dp MX 
xX 34, — = 
- v( vx = 2) at (° Yr) + 1 (= * ) 
= ¥ a 
dys DX 
2 2 df HX | 2 Wo? 
™ ri (srr) + AME r2 dz dz? at ii 
zr 
where 
ae 2 10d 
Beant jar e 
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and 


1) 
Ss iw — -+ —. 
row 


(c) Conditions for superposability of two motions.—Taking the two 
motions given by 


ea 1 dbs\- 
=(--")i, + Xt, + ( — I, 
“ fed si + ( r ar ie 
where s = 1, 2 corresponds to two motions 


Using the conditions of 
integrability (1.8) and (1.9), we obtain the following conditions for super- 
posability of two motions 


~ es ee +r) = r ia — 2y, es — 2X. 


Re 


= be}? (ter FEM) 42 5 (ter AEM) 
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(1.10) 
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dy dX, dy dX, dp, dX, 
+2){ dz dz 42) ow w J +2){ ow ow 
z er s Ff zr 
4.2 dp DX, 2 dy d°X, | 2 dy, °K, | 2 
row woz r Ww WIZ r oz dz r 
(1.11) 


We can easily show that the conditions (1.10) and (1.11) are both 
necessary and sufficient for the two axially symmetric flows of a non-New- 
tonian fluid corresponding to the two stream functions ¥%, and %, to be super- 
posable. 


2. CONDITIONS FOR SELF-SUPERPOSABILITY OF A FLOW: 


Using the conditions (1.10) and (1.11), it follows that the flow with 
stream function % is self-superposable if and only if, 


a(t) — a 


= 2v¢ |? J (* — — (* rev) 2d d ” 
shill i a" r dz 


ey wX X 
~pitey| +s r*J (Rr) 25 (3 ] 
. e 4 a 


“Re 7G (2.1) 
and 
X\ . dys dX db WX 
ae e[s(Wr) am Mz dz J+ 2) ( or ow 
‘ ae z ¢ s 4 
2 db 9°X 2 dob D°X 
- , dz + : |. (2.2) 
or woz raz 


3. PARTICULAR CASES 


(a) Newtonian fluids.—If in the conditions of superposability and self- 
superposability by putting ». = 0, we obtain the conditions of superposability 
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and self-superposability for a Newtonian fluid as obtained by Lakshmana 
Rao.* 

(b) Velocity vector lies in the meridian plane-——We can find the condi- 


tions of superposability and self-superposability in steady axi-symmetrical 
flows when the velocity vector g lies in the meridian plane, 


ie.y 


q = Url + Ugiz (3.1) 
or in terms of the Stokes’ stream function 
1d 1d 
dit (= i, +(-1™)iz (3.2) 


and the vorticity is given by 
= + EM. | (3.3) 


In this case the equations of motion are given by (1.14), (1.24) and 
(1.3a) with X = 0. 


The condition of integrability (1.8) becomes 
p rE 
rJ ( » ) 


= vEt} + 2v¢ i(? vm. “ (™ dy — 


Zz r 
1 dy (1 2 2m » 
~ a B¥— (584) + Fs (3.4) 


while the condition (1.9) is automatically satisfied. 
The condition of superposability (1.10) for two axi-symmetrical flows 
with stream functions %, and %, is given by, 


5 (%4 PEMA) 4.5 (2 EM) 
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FA dzZ r2 dz rz ar\ 7? 
2 dy, 2 (E* = 
Towzw vl. (3.5) 


The condition (1.11) is i SR satisfied. The condition for 
self-superposability of this flow is given by: 


ew 


aperm-a(tn) 


— (Bey +22 (EF \. (3.6) 


rodzor 


From the condition (3.5), we see that a rotational flow with suffix 2, is super- 
posable on an irrotational flow with suffix 1, if 


J (% r — 


22 fis(E) (9) 


Zz r 


La pay, 4 2 2 ae (3.7) 


~ 72 oz rozor\ r? 


Using the conditions (3.5) and (3.6), we can prove the following theorems. 


THEOREM 1.—Al/l irrotational motions are self-superposable and any 
two irrotational motions are superposable on each other. 

THEOREM 2.—Tne two flows will be superposable if the vorticity of each 
flow is proportional to rt. Also an axi-symmetrical flow will be self-super- 
posable if its vorticity is proportional to r. 

THEOREM 3.—If w,/r is constant along the stream lines of flow 2, while 
w/t is constant along the stream lines of flow 1, then the two flows are 
superposable if in addition the velocity components and vorticities satisfy the 
condition 

2 (UyY, + UVa) + F (Var, + Vw.) = 0. 
(where w, and a, are the vorticities of the two flows). 
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THEOREM 4.—An axi-symmetrical flow is self-superposable, if the vorti- 
city w is of the form w =f()r and if the stream function satisfies the 
condition 


fil 
EY => (5 EY). 
While the condition of integrability reduces to 


EN = “E> EN), 


r3 
Also a rotational flow with stream function ¥, is superposable on an 
irrotational flow. with stream function ¥,, if in general 
ws =f(p)r 
and 


EM, = 2r 2 (4 EM). 


The corresponding theorems for a Newtonian fluid have been given 
by Bhatnagar and Verma.* 


I am deeply grateful to Prof. P. L. Bhatnagar for suggesting the problem 
and for his kind guidance throughout this work. 
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INTRODUCTION 


THE kinetics of bromine addition to a-8 unsaturated acids under various 
conditions have been reported from time to time.!° The reaction in acetic 
acid is homogeneous, while in the non-polar solvent carbon tetrachloride 
is heterogeneous.2, A_ termolecular addition with increasing ‘bimolecular 
rate constants’ has been reported for allyl chloride and allyl bromide by 
Robertson et a/.'*% in acetic acid below 30°C. and the influence of electro- 
lytes on the addition reaction has also been studied. Nozaki and Ogg!® 
studied the reaction in acetic acid, while Juvala® used carbon tetrachloride 
where the normal electrophilic addition (allyl bromide faster than allyl 
chloride) has been reported. But Robertson etal. were unable to confirm 
Juvala’s observation even in an atmosphere of dry nitrogen. It was expected 
that a study of the thermal addition of bromine to allyl compounds—allyl 
chloride, allyl bromide, 2-methyl allyl chloride and allyl propionate—in 
acetic acid and in acetic acid-carbon tetrachloride solvent mixtures. at 
different temperatures would lead to a better understanding of both struc- 
tural and solvent influences. 


EXPERIMENTAL 


Materials Used.—Bromine, acetic acid and carbon tetrachloride were 
purified as in earlier communications.’® Allyl chloride and allyl bromide 
were purified by the method of Vogel.4® Calcium bromide was used for 
drying allyl bromide, while 2-methyl allyl chloride was purified by the same 
method as allyl chloride. Allyl propionate was washed with a saturated 
solution of sodium bicarbonate and then with water, and dried over a 
prolonged period by calcium chloride and fractionated. Periodic purifica- 
tion before use was found to be necessary with all allyl halides. The com- 
pounds were stored in amber-coloured bottles in the dark, in desiccators 
over fused anhydrous calcium chloride. 
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Allyl chloride .. 45°C. (758 mm.). 

Allyl bromide .. 70°C. (759 mm.) 

Allyl propionate .. 123°C. (760 mm.) 

2-methyl allyl chloride .. 722°C. (757 mm.) 

Acetic acid .. 118°C. (760 mm.) [n]p*° = 1-3672 
Carbon tetrachloride ..  76°8° C. (760 mm.) [n]p*° = 1-4539 


Solvent mixtures were prepared by adding together appropriate volumes 
of each component, taking the usual precautions. These mixtures were 
used for preparing the standard solutions of olefine and bromine. All 
weighings were carried out only in weak diffused light, due to the sensitivity 
of allyl compounds to light. Typical experimental results are presented 
in Figs. 1, 2, 3 and 4. 
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Fics. 1-4. Fig. 1. Allyl Bromide. 89%HOAc. Fig. 2. Allyl Chloride (M/80). 
Fig. 3. Allyl Propionate (M/160). Fig. 4. 2-M-thyl Allyl Chloride (M/80). 


In all cases, the reaction is homogeneous, differing from the olefine- 
acids, which in pure carbon tetrachloride were found to show heterogeneous 
addition. Since Robertson et al. *° had reported a termolecular addtion, the 
termolecular rate constants were calculated by the equation 


Ke; (ap ah 
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DISCUSSION OF RESULTS 


Robertson has reported that the order of the reaction between allyl 
chloride and bromine, at concentration 0:04 to 0:02 M, was 3-2, some- 
what higher than the values 3-0 to 3-1 obtained for other compounds over 
this concentration range, and this higher value has been attributed possibly 
to autocatalysis. He considered that Nozaki and Ogg’s data show an order 
3-3 and 3-7. The present results too appear to indicate a third order, the 
order with respect to olefine being invariably one. The increasing bimole- 
cular constants obtained by them with the progress of the reaction clearly 
indicated autocatalysis. Retardation by the dibromides may not be pro- 
nounced in acetic acid, as is shown by our investigation. The hydrogen 
bromide formed acts as a catalyst in the case of olefine acids, but the effect 
may not be as pronounced for allyl compounds as for acids, since substi- 
tuents may be less deactivating in the former group of compounds. An 
order greater than two is consistent with decreasing K, values, while in- 
creasing K, values may indicate an order less than three in the absence of 
disturbing influence on the reaction, as autocatalysis. 


We find that in all cases, bimolecular constants decreased with the 
progress of the reaction up to about 40° completion of the reaction 
but fairly good termolecular constants have been obtained. 


TABLE I 


Bimolecular constants (K,) and termolecular constants (K;) for the reaction 
between bromine and allyl chloride, and bromine ard allyl brcmide 
in the solvent mixtures (Acetic acid-Carbon tetrachloride) at 28° C. 


(Concentration of reactants = M/80) 





Value for % completion of the reaction 
Olefine Solvent % Ky K; 
10 40 10 40 











Allyl chloride -- 100 HOAc 1-347 1-088 454-9 453-4 
~ - w@ 0-8081 0°6311 273-0 269-3 
~~ a « 0-4136 03194 139+7 136°3 
« SS ww 0-1022 00-0833 34°51 35-56 

Allyl bromide co TP a 1-1117 =0-9155 =. 395-1 390°4 
‘ — « 0-5926 00-4618 200-2 192-5 

80 ly, 0-2778  0-2265 93-83 96-63 
@ « 0-0773 0-0591 25-66 25-23 





In all cases Kg values are presented as moles per litre per minute and K, values as (moles per 
litre)? per minute. 
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The tendency for the termolecular constants to increase is also noted 
in the later stages of the reaction at lower concentration of allyl propionate 
and 2 methyl allyl chloride. This can be ascribed to catalysis by hyd-ogen 
bromide formed during the course of the reaction and also the gradual 
dominance of the bimolecular reaction over the termolecular one. Further, 
it is seen that the termolecular constants appear to show an inverse variation 
with concentration, indicating clearly a composite mechanism of bimolecular 
and termolecular types. The bimolecular reaction may not be so pronounced 
during the initial stages at the concentration investigated. 


TABLE II 


Effect of concentration changes on Ky values (Kj X 10~*) 


(Solvent: Acetic acid-Carbon tetracl.loride mixture. Temperature 28°C.) 





Concentration 








Olefine Solvent % 
M/80 M/160 M320 
2-methyl allyl chloride 100 HOAc +“ oe 218-30* 
~ 80 os 43-32 48-04* 
™ @ ws 9-871 10-87 12-49 
™ 40 ww 2-45 2+425 2°734 
Allyl propionate @ 17-29 23-04 
” eS « 7°285 8-668 
‘i . _ « 2205 2°779 





* At 10% completion of the reaction. 


The reaction shows a pronounced sensitivity to added water (Fig. 5). 
In 99% acetic acid (1% water) K, increases initially and with increased 
amounts of water K, tends to remain constant, whereas in pure 
acetic acid, K. decreases. This may be the resultant of two changes: a 
change in the polarity of the environment and formation of the hydrogen 
bromide by interaction of water with bromine. The termolecular reaction 
thus becomes a bimolecular one and the increasing K, values arise from 
autocatalysis, 
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An analysis of earlier work clearly indicates that in the addition of 
bromine, increased dilution of reactants, rise of temperature, or small 
increase in water content of the solvent results in a gradual dominance of 
a second order rate with a higher heat of activation. The tendency of K, 
values to increase with the progress of the reaction is greater in the more 
dilute solutions and also at the higher temperature. We find the same trend. 


P20 
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Fig. 5. Allyl Chloride (M/80). 


The change of order with dilution may now be analysed. When the 
halogens and the olefines react, the first result is the apparent formation of 
a 1:1 molecular complex. The termolecular process may then be written 
thus : 


K, Br, _ 
2 Ks solvent 
C (hydro- 
—+ ABr.Br+C —— 
K. eating 


The rate of formation of the complex (A, Br,) will be proportional to the 
encounters between olefine and bromine, and the rate of disappearance -will 
be the sum of two factors: the decomposition of (A, Br.) into A and Br, 
and their disappearance in the actual chemical transformation. 


—_— = K, (A) (Br,) — K, [A, Br,] — Kg [ABr,] [C] 


Since for the steady state d[A, Br,]/dt = 0, the rate of reaction is equal to 


K,-K,- (A) Br.) (C) 
es eo, le 


K, [ABr,] [C] = 














Kinetics of the Olefine Bromine Reaction—12 169. 


At low concentrations of bromine.(C = Br,), the reaction rate is thus pro- 
portional to (A) (Br,)* and at high concentration, the rate is proportional 


to (A) (Bry). 


According to Robertson’ in the concentration region (M/20) — (M/80) 
in acetic acid, the reaction is termolecular. Even at lower concentration 
in acetic acid, carbon tetrachloride solvent mixtures, the reaction appears 
to be termolecular. The rate thus appears to be proportional to (A) (Br,)?. 
The decrease in order with dilution in acetic acid and the apparent bimole- 
cular constants as well as the observed bimolecular reaction in solvents more 
polar than acetic acid has then to be accounted for. With increasing dilution 
the place ‘of a second bromine molecule is apparently taken by the acetic 
acid molecules and in the more polar solvents this is even more so. Hence, 
the reaction 


[ABr.} — ABr* + Br 
| HOAc 
ABrOAc + H* 
ABr+ + Br — ABrBr. 


The primary association of olefine and bromine thus becomes the rate deter- 
mining step. The bromine addition in methyl alcohol is bimolecular*® 
because of the higher dielectric constant which facilitates polarisation of 
the bromine molecule and the greater reactivity of (OMe)- compared to 
the (OAc)-. Catalysis by HBr may not be so pronounced for a termolecular 
reaction. It has been observed by Nozaki and Ogg that more bromide 
ions are formed at lower concentration of bromine in acetic acid during 
the reaction. According to Andrews® in aromatic halogenation in hydro- 
xylic solvents, the place of bromine is replaced by the solvent. Hence ‘in 
the equation: 


= = K, (A) (Br) + Ks (A) (Br,)? 


K, and K, correspond to the functioning of acetic acid and bromine as 
electrophilic reagent in promoting the anionisation of Br, from the complex 
(A, Bry). 


Effect of temperature, substituents and solvents can be assessed from 
Tables III and IV. It can be seen from Table III that a plot of log K against 
1/T will not be a straight line indicating a mechanistic change with temperature; 
tle energy of activation increases with temperature. 
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TABLE III 


Effect of temperature : 2-methyl allyl chloride and bromine reaction 
Concentration of reactants: M/1€0 


(Solvent: Acetic acid-Carbon tetrachloride mixture) 











So'vent Tempe- K, E Average PZ 
y 4 rature K.cal. E 
80 HOAc ae 28 4332 . 
1-82 
35 4642 2°58 5+1x 108 
} 3-34 
40 5041 
60 HOAc i 28 1081 — 
ae | 
35 1145 1-87 4-0x 10? 
2°37 
4G 1218 
40 HOAc es 28 242-5 
0-49 
35 247-0 0-83 3-0x10 
1-17 
40 254-2 
TABLE IV 


Effect of substituents 
(Solvent: 60% HOAc, 40% CCI) 








Allyl Allyl 2-methyl Ally! 
chloride bromide allyl chloride propionate 
M/80 M/80 M/160 M/160 
E. K. Cal. (Average) - 1-08 3°82 1°87 4°69 
Log,,PZ . 2 2-38 2-60 4-54 





Changes in the Arrehenius parameters follow the same trend except 
for allyl bromide, where there seems to be a decrease in the apparent order 
of the reaction with increase in carbon tetrachloride content of the medium. 
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The values of the frequency factor decrease in all cases with more of carbon 
tetrachloride but the decrease is not so great for allyl bromide. The change 
in the value of E and PZ can arise from an increase participation of a ter- 
molecular mechanism, with increasing carbon tetrachloride content which 
in its turn decreases the apparent polarity of the solvent. A study of absorp- 
tion spectral! of bromine in acetic acid and carbon tetrachloride indicates 
clearly the solvation effects in hydroxylic solvents. The acetic acid may 
be forming a semi-rigid cage around the bromine molecule with the resultant 
electrostatic polarisation of the halogen molecule. The extent of this 
decreases with the decreasing polarity of the solvent. This view is supported 
by the linear plots obtained for log K; against (D—1)/(2D+1) for this 
dipole reaction (Fig. 6). 


Using the conventional rate equation, we have K, = PZ,e-™/8T 
where Z, is the collision rate at unit concentration. For a termolecular 
reaction Zs, is of the order 10° to 10!%, 1? moles~® sec.—! and the values of P 
is of the order 10-5 to 10-'° in reaction involving polyatomic molecules,®,’ 
low values of the frequency factor have to be expected and the results bear 
this out. 


Migration of groups as well as the relative stability of carbonium ion 
which in its turn is affected by the presence of substituents cannot be 
excluded in the addition of electrophilic reagent such as HOC]! to the allyl 
halides.©5 Now the bromine addition to allyl halide may be pictured thus: 


+ 8- + 
BrCH,CHCH,x —> CH,;——-————CH——CH,x 
P 
‘ ? 
ra 
Np’ &- 
11 
x & «x 8 
i % ff 
YQ Pa 
*, + ri 








BrCHy—CH———CH, = s'cu;——cH——cu* 
a ¢ 


When X = Br, the possibility of the chances of such a structure will be 
increased when compared with allyl chloride, i.e., such complete equivalence 
may not be expected for allyl chloride as for allyl bromide for the bromine 
addition reaction. This can account for the difference between allyl chloride 
and allyl bromide, in the latter compound the decrease in the frequency 
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factor is less. Two products of bromine addition to allyl chloride aré 
possible, namely, CH,BiCHBrCH,Cl and CH,BrCHCICH,Br while for 
allyl bromide the sole product of the bromine addition is CH,BrCHBrCH,Br 
assuming only the dibromide formation. This aspect can only be confirmed 
with tracer techniques. Further, the halogen hyperconjugation will be 
greater for allyl chloride than in allyl bromide. Electron releasing groups 
in the 2-position should stabilise the carbonium ion and allow it to last 
for a longer life before reaction with the nucleophilic part. In this way, 
the reaction through the symmetrical intermediate may be favoured for 
2-methyl allyl chloride and as a result the products may be CH,BrC (CH,) 
CICH.Br and CH,BrC (CH;) BrCH.Cl. Complete equivalence of the 
halogen may not be expected since chlorine is more electronegative than 
bromine and the products formed will not be in equal proportions. This 
is considered in a later communication. 
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Fic. 6. Allyl Bromide (M/160). 


In the series, 2-methyl allyl chloride, allyl propionate, allyl chloride 
and allyl bromide, in pure acetic acid, the reactivity decreases in the order, 
but the addition of carbon tetrachloride brings about a change, allyl pro- 
pionate becoming faster than 2-methy] allyl chloride, as is seen from Table V. 


The relative reactivity of allyl chloride and bromide has been attributed 
by Robertson etal. to the two modes of reaction: (1) due to hydrogen 
hyperconjugation and (2) due to the inductive effect of the halomethyl 
group. This may be extended, in the present case, to 2-methyl allyl chloride. 
The methyl group facilitates the initiation of the addition at the ethenoid 
carbon atom remote from the methyl group (Carbon — 1) +I effect of 
methyl group thus augments the C—H_ hyperconjugation of the methyl 
group as well as that of the chloromethyl group which in its turn is effected 
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TABLE V 


Influence of substituents and solvent composition (concentration M/160) 
(Solvent: Acetic acid-Carbon tetrachloride mixture) 





~=—_—__ 





Solvent Allyl Allyl 2-methyl Allyl 
% chloride bromide allyl chloride propionate 

80 HOAc ‘ 149-5 110-1 4332 2304 
@ « fe 38-58 30-05 1087 866-8 
® » rr 6 242°5 277°9 





by —I effect of the chlorine, and the —I effect can initiate the addition at 
carbon —2 as shown below. 


f\ 
H,C—H 
CH 
L\ f/¥41 ' 
CH gam CH game 
1 2 1 2 
CH, 
CHCl «f+ 
AN Cl 
H 


Further, the polarising forces of the methyl and chloro-methyl groups are 
—0:39 and +0-36C.G.S. units!* and this indicates that initiation of addi- 


tion at carbon (1) is more important. 


The rate of addition of bromine to propylene is greater than 
allyl chloride by a factor greater than 10%. By reference to the reactivities 
of ethylene and allyl chloride, the reactivity of bischloro-ethylene is given 
as 1:0°019 :0-00024. 2-methyl allyl chloride is found to be faster than 
allyl chloride by a factor of 30. Hence, 2-methyl allyl chloride adds bromine 
considerably rapidly when compared with bischloro-ethylene, the rate factor 
being greater than 10%, as must be expected due to the powerful electron 
withdrawing influence of the halogen substituent. Hence the slow addition 
of bischloro-methyl ethylene cannot be attributed to steric hindrance, even 
if it operates, it will be of minor importance. 


In allyl propionate —I effect of ‘CO’ may not have an appreciable- 
effect.on the ethylenic system for it is far removed from the ethenoid centre, . 
Hence, allyl propionate must be faster than allyl halides. Further, it gives 
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a clear evidence for a composite mechanism—both bimolecular and ter. 
molecular—are operative and the bimolecular reaction slowly changes over 
to a termolecular one with a decrease in the polarity of the solvents. This, 
in turn, is reflected in the energy of activation and frequency factor. Unlike 
a-8 unsaturated acids, addition to allyl compounds does not indicate 
any high degree of autocatalysis. 


Hence, it is clear that in comparing the influence of substituents, the 
velocity constant alone is inadequate. Further, if one compares the reacti- 
vity using the velocity constant, there is the implicit assumption that fre- 
quency factor is identical for all the olefines and this assumption can, by 
no means, be justified. Substituents may not influence equally in the gas 
phase and in solution and may also vary with the solvent (cf. 1). 


SUMMARY 


The thermal addition of bromine to allyl compounds has been studied 
in acetic acid-carbon tetrachloride solvent mixtures at various temperatures. 
The influence of solvents and substituents on the reaction mechanism has 
been discussed. 
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